Limits on electron quality in suspended graphene due to flexural phonons 
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The temperature dependence of the mobility in suspended graphene samples is investigated. In 
clean samples, flexural phonons become the leading scattering mechanism at temperature T > 10 K, 
and the resistivity increases quadratically with T. Flexural phonons limit the intrinsic mobility down 
to a few m 2 /Vs at room T. Their effect can be eliminated by applying strain or placing graphene 
on a substrate. 



Introduction. — The properties of isolated graphene 
continue to attract enormous interest due to both its ex- 
otic electronic properties [l[ and realistic prospects of 
various applications 0]. It has been found that the in- 
trinsic mobility fx of charge carriers in graphene can ex- 



ceed 20 m 2 /Vs at room temperature T Jjj, 0], which is 
the absolute record. So far, such high values have not 
been achieved experimentally, because extrinsic scatter- 
ers limit it. The highest \i was reported in suspended 
devices [1, Q and could reach ~ 12 m 2 /Vs at 240 K Q. 
This however disagrees with the data of Ref. [f| where 
similar samples exhibited room T /x close to ~ 1 m 2 /Vs, 
the value that is routinely achievable for graphene on a 
substrate. 

In this Letter, we show that flexural phonons (FP) 
are an important scattering mechanism in suspended 
graphene and the likely origin of the above disagreement, 
and their contribution should be suppressed to allow 
ultra high fi. Generally, electron-phonon scattering in 
graphene is expected to be weak due to very high phonon 
frequencies Q. However, in suspended thin membranes, 
out of plane vibrations lead to a new class of low energy 
phonons, the flexural branch In an ideal flat sus- 

pended membrane symmetry arguments show that elec- 
trons can only be scattered by two FP simultaneously 
0, El ■ As a result the resistivity due to FP rises rapidly 
at high T where it can be described as elastic scattering 
by thermally excited intrinsic ripples 12j. 



We analyze here the contribution of FP to the resis- 
tivity, and present experimental results which strongly 
support the suggestion that FP are a major source of 
electron scattering in suspended graphene. This intrin- 
sic limitation to the achievable conductivity of graphene 
at room T can be relaxed by applying tension, which 
modifies both the phonons and their coupling to charge 
carriers. 

Model. — Graphene is a two dimensional membrane, 
whose elastic properties are well described by the free 



energy 
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where K is the bending rigidity, A and /x are Lame co- 
efficients, h is the displacement in the out of plane di- 
rection, and Uij = 1/2 [d{Uj + djUi + (dih)(djh)] is the 
strain tensor. Summation over indices in Eq. CP is 
implied. Typical parameters for graphene [13l4l5| are 
K « 1 cV, and /x w 3 A « 9 eV A -2 . The density is 
p = 7.6 x 10~ 7 Kg/m 2 . The velocities of the longitudinal 
and transverse phonons obtained from Eq. ([1]) are = 

~ 2.1 x 10 4 m/s and v T = ~ 1.4 x 10 4 m/s. 

p ' V p 

The FP show the dispersion 



a q' 



(2) 



with a 



4.6 x 10~ 7 m 2 /s. 



Suspended graphene can be under tension, either due 
to the electrostatic force arising from the gate, or as a 
result of microfabrication. Let us assume that there are 
slowly varying in plane stresses, Uij(r), which change lit- 
tle on the scale of the Fermi wavelength, kp 1 , which is 
the relevant length for the calculation of the carrier re- 
sistivity. Then, the dispersion in Eq. ([2]) is changed into: 



«f (?) = |q|J-|q| 2 + ^(r) + (3) 

q ]j p p p |qr 

The dispersion becomes anisotropic. For small wavevec- 
tors, the dispersion is linear, with a velocity which scales 
as y/u, where u is strain. 

The coupling between electrons and long wavelength 
phonons can be written in terms of the strain tensor. 
On symmetry grounds, we can define a scalar potential 
and a vector potential which change the effective Dirac 
equation which describes the electronic states 



V(r) = go [u xx (f) + u yy (f}] 

A(f) = ^ [u xx (r) - u yy {f)},u xy (r) 



(4) 
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where go ~ 20 — 30 eV is the bare deformation potential 
fl(| . a 1.4 A is the distance between nearest carbon 
atoms, P = -d\og(t)/dlog(a) f« 2 - 3 [H, and t w 
3 eV is the hopping between electrons in nearest carbon 
7r orbitals. 

Linearizing Eq. (TTJ) and expressing the atomic displace- 
ments in terms of phonon creation and destruction op- 
erators, and using Eq. (j4]) and the Dirac Hamiltonian 



for graphene jl[ we can write the full expressions for the 
coupling of charge carriers to longitudinal, transverse and 
FP, without and with preexisting strains. 

Calculation of the resistivity. — We assume that the 
phonon energies are much less than the Fermi energy, 
so that the electron is scattered between states at the 
Fermi surface. After some algebra, the scattering rate 
due to FP, including a constant strain, u, is 
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FIG. 1: (color online), a) Two phonon diagram which de- 
scribes electron scattering, b) Kinematics of the process, and 
variables used in Eq. ©■ The circle denotes the Fermi sur- 
face. 



where vp is the Fermi velocity and kp the Fermi momcn- 
turn, [D{K)f = [g(K)] 2 [l-K 2 /(2kp) 2 ] + {f5hvp) 2 /(Aa 2 ) 
is the generalized deformation potential, including the 
contribution of the screened scalar potential g{K) 
go/e(K) and gauge potential, uj q 



uv L q ,L is 



the phonon dispersion, Eq. ([2]) in the isotropic approx 
imation, and n q is the Bosc-Einstein distribution func- 
tion. The diagram described in this calculation, and 
the variables K = |K|, q = |q| and Q = |Q| are 
shown in Fig. [TJ The static dielectric function is e(K) = 
1 + e 2 N(k F )/(2e K), where N(k F ) = (2k F )/ (irhvp) is 
the density of states. At kp the screened scalar potential 
g = go/£(kp) « go/8 ~ 3eV is in good agreement with 
ah initio calculations (20| . 

The relevant phonons which contribute to the resis- 
tivity are those of momenta |q| > 2kp. This scale 
allows us to define the Bloch-Griineisen temperature, 
ksTBG = fiaj2k F - Neglecting first the effect of strain, 
we find: 
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respectively for in-plane longitudinal (L) and transverse 
(T) and for FP (F 1 ), where the temperature is in Kelvin 
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and the electron density n is expressed in 
Close to room T we arc in the regime T T^q' 1 ,r for 
all concentrations of interest. The corresponding tem- 
perature for FP in the presence of a uniaxial strain, u is 
Tbg — 2%^/u^n K. Our focus here is on the experimentally 
relevant high-T regime. 

In systems with strain, the phonon dispersion relation, 
Eq. (p|). shows a crossover between a regime dominated 
by the to another where the strain becomes irrelevant, at 
q* w vpsfXij a. The range of integration over the phonon 
momenta in Eq. ([5]) is limited by hu> q < fc^T, and kp < q. 
In addition the theory has a natural infrared cutoff with a 
characteristic momentum q c below which the anharmonic 
effects become important |2l|. Defining qp as huj qT = 



ksT, the scattering rate in Eq. shows three regimes 
in which (i) strain is irrelevant and max(q*,q c ) <C kp, 
(ii) strain is small and relevant phonons combine linear 
and quadratic spectrum for rn&x(kp,q c ) < q* < qp, (in) 
strain is high and determines the scattering rate for qp <C 
q*. We finally obtain: 
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where D 2 = g 2 /2+(Phvp) 2 /(4a 2 ), and the infrared cutoff 
Huj c is related to max(q*, q c ). For comparison we give also 
the contribution from in-plane phonons, 
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The T dependence of the scattering due to FP is more 
pronounced than that due to in-plane phonons, and it 
dominates at high enough T. In the limit of irrelevant 
strains, max(q*,q c ) <C kp, the crossover temperature is 



tl,t{ t *\ 
r F (T*) 



= 1 T*(K) « 57 x n(10 12 cm- 2 ). (9) 
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FIG. 2: (Color online), (a) Contribution to the resistiv- 
ity from flexural phonons (blue full line) and from in-plane 
phonons (red dashed line), (b) Resistivity for different strain. 
The in-plane contribution (broken red line) shows a crossover 
from a low to a high-T regime. In both cases, the electronic 
concentration is n ■ 



When T* < Tgg this crossover does not occur and scat- 
tering by FP dominates also at low temperatures. At 
finite strain max(fcj?,q c ) <C q* <C qr we obtain 



tl,t{ t *) 
tf{T*) 



1 T*(K) ss 10 6 m. 



(10) 



In the absence of strains, the crossover shown in Eq. © 
implies that the room T mobility is limited by FP for den- 
sities below 10 13 cm -2 . Strains reduce significantly the 
effect of FP, so that, in the presence of strain, the mobil- 
ity is determined by the scattering by in-plane phonons, 
see Eq. (fTOj) . 

The contribution to the resistivity from the different 
phonon modes can be written, using the expressions for 
the scattering rate as 



Pi{n,T, u) 



e 2 v 2 F N(k F )Ti(n 7 T, u) 
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where the index i label the phonon mode. Results for the 
resistivity in different regimes are shown in Fig. [5] 

The results of Eqs. and ([TTj) can be extended to 
bilayer graphcnc. The main differences arc: (i) the kine- 
matics of the two phonon scattering are the same as in 
Fig. [TJ except that the overlap between the electronic 
states |k) and |k+K) is modified; (ii) the density of states 
is constant N(e) = N = t±/(nh 2 v' F ), and the screened 
scalar potential g is replaced by (2eogokp)/[e 2 N(eF)] = 
(2ireogokFfi 2 v F ) / (e 2 t±), where t± is the hopping between 
layers, and e is the electric charge; (iii) the Fermi velocity, 
which determines the coupling to the gauge potential is 
Ivpftkp /t±. The Fermi velocity and the density of states 
also change the expression for the resistivity, Eq. (jTTJ) . 

Experimental results. — We have fabricated two- 
terminal suspended devices following the procedures in- 
troduced in Refs. [H, @. Typical changes in the resis- 
tance R as a function of the gate induced concentration 
n are shown in Fig. 3a. The as-fabricated devices ex- 
hibited jU ~ 1 m 2 /Vs but, after their in situ annealing 



by electric current, fj, could reach above 100 m 2 /Vs at 
low T. To find fx, we have used the standard expression 
R = R + {I / w){\ I nen) where R describes the contact 
resistance plus the effect of neutral scatterers, and both 
Rq and fi are assumed n- independent [H, H| ■ Our devices 
had the length I ps 1 — 2 /jm and the channel width w 
of 2 - 4 fjm (see the inset in Fig. 3b). At T > 100 K, 
the above expression describes well the functional form 
of the experimental curves, yielding a constant (i over 
the whole range of accessible n, if we allow Ro to be 
different for electrons and holes. This is expected be- 
cause of an n — p barrier that appears in the regime of 
electron doping due to our p-doping contacts At 
T < 100 K, the range of n over which the expression fits 
the data rapidly narrows. Below 20 K, we can use it only 
for n < ±10 10 cm -2 because at higher n we enter into 
the ballistic regime (the mean free path, proportional to 
/in 1 / 2 , becomes comparable to I). In the ballistic regime, 
graphene's conductivity a is no longer proportional to 
n |,|] and the use of [i as a transport parameter has 
no sense. To make sure that /i extracted over the nar- 
row range of n is also correct, we have crosschecked the 
found /i against quantum mobilities inferred from the on- 
set of Shubnikov-de Haas oscillations 0, Q . For all our 
devices with \x ranging from ~ 1 — 100 m 2 /Vs, we find 
good agreement between transport and quantum mobili- 
ties at liquid-helium T. Fig. 3b shows the T dependence 
of fj.. It is well described by the quadratic dependence 
1/jtx = 1 /fi(T -> 0) + 7T 2 . Surprisingly, we find the coef- 
ficient 7 to vary by a factor of ~ 2 for different devices, 
which is unexpected for an intrinsic phonon contribution. 
Such variations are however expected if strain modifies 
electron-phonon scattering as discussed below. Note that 
H falls down to 4 - 7 m 2 /Vs at 200 K (sec Fig. 3b) and 
the extrapolation to room T yields \i of only 2 — 3 m 2 /Vs, 
which is significantly lower than the values reported in 
Ref. but in agreement with Ref. 

Discussion. — The density independent /1 indicates 
that experiments are in the non-strained regime, where 
1/rp - T 2 /k F and p F - T 2 /n. Here FP completely 
dominate and the coefficient 7 defined above is given 

by 7 w r , D k B -> i hi(4 a ^), where the infrared cutoff is 
the only free parameter [22| . Experiment gives 7 ~ 
6.19 x 10~ 6 Vs/ (mK) 2 for the sample with lower mobility 
and 7 ps 3.32 x 10~ 6 Vs/(mK) 2 for the higher mobility 
one. Neglecting the logarithmic correction of order unity, 
the analytic expression gives 7 ps 3x10~ 6 Vs/(mK) 2 
without adjustable parameters. 

The difference between the two samples may be un- 
derstood as due to a different cutoff under the logarithm 
due to strain. In non-strained samples there is a natural 
momentum cutoff q c ps 0.1 A below which the har- 



monic approximation breaks down [2l| . Strain increases 
the validity of harmonic approximation, making q c strain 
dependent, thus explaining different cutoff at different 
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FIG. 3: (color online), (a) Electron transport in suspended 
graphene. Graphene resistivity q = R(w/l) as a function of 
gate-induced concentration n for T = 5, 10, 25, 50, 100, 150 
and 200 K. (b) Examples of fi(T). The T range was limited 
by broadening of the peak beyond the accessible range of n. 
The inset shows a scanning electron micrograph of one of 
our suspended device. The darker nearly vertical stripe is 
graphene suspended below Au contacts. The scale is given by 
graphene width of about 1 fim for this particular device. 

strain. A rough estimate of the expected strains is ob- 
tained by comparing q c ps 0.1 A with q* = VL\/H/a, 
which gives u ~ 10~ 4 - 1CT 3 . Such small strain can be 
present even in slacked samples (where strain induced by 
gate and T is negligible) due to, for example, the initial 
strain induced by the substrate and remaining unrelaxed 
under and near metal contacts. A complete theory would 
require the treatment of anharmonic effects, which is be- 
yond the scope of the present work. The data in Q show 
higher mobilities than those in Fig. [3] A fit to this data 
using Eq. suggests the sample being under strain. 

Conclusions. — The experimental and theoretical re- 
sults presented here suggest that FP are the main mech- 
anism which limits the resistivity in suspended graphene 
samples, at temperatures above 10 K. Scattering by FP 
involves two modes, leading to a T 2 dependence at high 
temperatures, with mobility independent of carrier con- 
centration. These results agree qualitatively with classi- 
cal theory assuming elastic scattering by static thermally 
excited ripples [l2|. Quantitatively, one of our main re- 
sults is that in devices with negligible strain the mobility 
does not exceed values of the order of 1 m 2 V _1 s -1 at 
room T, that is, FP restrict the electron mobility to val- 
ues typical for exfoliated graphene on a substrate. 

The dispersion of FP changes from quadratic to linear 
when the sample is under tension. As a result, the in- 
fluence of FP on the transport properties is suppressed. 
The T dependence of the mobility remains quadratic, but 
it decreases linearly with the carrier concentration. Im- 
portantly, applying rather weak strains may be enough 
to increase dramatically the mobility in freely suspended 



samples at room T. 

A very recent theory work j23[ has also addressed the 
role of FP on electron transport. Insofar as the two anal- 
ysis partially overlap, the results are in agreement. 
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